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ABSTRACT 

We propose a simple model for the formation of dwarf spheroidal galaxies, in which stars are assumed to 
have formed from isothermal gas in hydrostatic equilibrium inside extended dark matter halos. After expelling 
the leftover gas, the stellar system undergoes a dynamical relaxation inside the dark matter halo. These models 
can adequately describe the observed properties of three (Draco, Sculptor, and Carina) out of four Galactic 
dwarf spheroidal satellites studied in this paper. We suggest that the fourth galaxy (Fornax), which cannot be 
fitted well with our model, is observed all the way to its tidal radius. Our best fitting models have virial masses 
of ~ 10 9 Mq, halo formation redshifts consistent with the age of oldest stars in these dwarfs, and shallow 
inner dark matter density profiles (with slope 7 ~ -0.5 . . .0). The inferred temperature of gas is ~ 10 4 K. In 
our model, the "extratidal" stars observed in the vicinity of some dwarf spheroidal galaxies are gravitationally 
bound to the galaxies and are a part of the extended stellar halos. The inferred virial masses make Galactic 
dwarf spheroidals massive enough to alleviate the "missing satellites" problem of ACDM cosmologies. 
Subject headings: dark matter — early universe — galaxies: dwarf — galaxies: formation — methods: A/-body 
simulations 



1. INTRODUCTION 

The popular A cold dark matter (ACDM) cosmological 
concordance model has been very successful in describing the 
properties of the universe on large and intermediate scales, but 
has encountered numerous difficulties on smaller, galactic and 
sub-galactic scales (e.g. Tasitsiomi 2003). Among the most 
obvious discrepancies is the prediction that there should be a 
factor of 50 more satellites orbiting in the Milky Way halo, 
with the masses comparable and larger than the conventional 
estimates for Galactic dwarf spheroidals (dSphs), than is ac- 
tually observed — th e so-called "missing satel lites" pr oblem 
(Moore et afl ll999l) . IStoehr et all J2002l> and lHavashi et"aTI 
(2003) gave a possible solution to this problem. They argued 
that the properties of Galactic dSphs are consistent with them 
being embedded in very extended dark matter (DM) halos, 
with total masses comparable to the masses of the largest sub- 
structure predicted to populate the Milky Way halo by ACDM 
cosmologies. 

There have also been cosmology-independent arguments 
in favor of much more massive dSphs than the conventional 
~ 10 7 M^ estimates ba se d on th e "mas s follows light" as- 
sumption (Mateol [T998l) . lLokasI ( 120021) modeled the line- 
of-sight velocity dispersion profiles for Fornax and Draco 
under the assumption that the anisotropy parameter j3 = 
l-(a,/a r ) 2 is a constant, and derived the total masses for 
these galaxies of ~ (1-4) x 10 9 M Q . (Here 07 and a r are 
the tangential and radial velocity dispersion, respectively.) 
From fitting a two pa rameter family of spherical models of 
Wilki nson et alJ (120021) to th e observed line-of-sig ht velocity 
dispersion profile in Draco, iKlevna et al.l ll2002l) concluded 
that the mass-to-light rat io in this dwarf increases outwards. 
lOdenkirchen et al J ( 120011) used large-field multicolor photom- 
etry of an area of 27 square degrees in the vicinity of Draco 
to show that this galaxy has very regular stellar isodensity 
contours down to their sensitivity limit (0.003 of the cen- 
tral surface brightness), with no signs of tidal interaction 
with the Milky Way. This implies that Draco has more DM 
in its outskirts than in its center and a practically uncon- 



strained up per limit on total mass. [B urkert & Ruiz-Lapuentd 
( 119971) and iMashchenko. Carignan. & Bouchard! ( 12004 pro- 
posed two different models (supernovae la heating of the ISM 
and radiation harassment by the ionizing photons from the 
host galaxy, respectively) to explain complex star formation 
history of some dSphs. The important ingredient of both 
scenarios is the requirement for dSphs to be very massive 
(> 10 9 MQ) so that the fully ionized ISM would not escape 
to intergalactic space. 

If it is indeed the case that many Galactic dSphs have mas- 
sive, 1 9 — 1 1 M Q DM halos (which would translate into 
their apparent tidal radii being > 5 angular degrees in the sky), 
then the observed distribution of stars in these galaxies should 
have been entirely shaped by internal processes, and not by 
the tidal field of the Milky Way. The current stellar density 
in Galactic dSphs is extremely low (they are DM dominated 
even at the center), with the estimated dynamical relaxation 
time for stars being much larger than the Hubble time (as- 
suming that the DM halo is sufficiently smooth at the center). 
It is reasonable to assume then that the distribution of stars 
in dSphs has stayed unchanged since the initial star formation 
epoch. 

lLakel lfl990) noted that the core line-of-sight velocity dis- 
persion in Draco and Ursa Minor (~ 10 km s -1 ) is of the same 
order as the sound speed in a diffuse cooling cloud of pri- 
mordial composition. He argued that a simple model of stars 
forming from isothermal gas in hydrostatic equilibrium in a 
DM halo with the Plummer density profile provides a reason- 
ably good description of the observed properties of these two 
dSphs. His conclusion was that the apparent cutoff radius in 
the distribution of stars in Draco and Ursa Minor has nothing 
to do with galactic tides. 

Here we propose a simple star formation model for dSphs. 
We assume that stars formed from isothermal self-gravitating 
gas in hydrostatic equilibrium inside an extended DM halo. 
We simulate the relaxation of the stellar cluster after the left- 
over gas is expelled using a n A/-body cod e. We significantly 
improve upon the model of lLa 

S ©90) by (1) considering 
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self-gravity of gas and stars, (2) allowing for dynamical relax- 
ation of stars after expelling the leftover gas, (3) introducing a 
lower density cutoff for star formation, and (4) using properly 
normalized ACDM DM halos. We show that our models are 
consistent with the existing observations on three out of four 
dSphs studied in this paper, if one assumes that these dwarfs 
have very massive (~ 10 9 M Q ) DM halos. 

The paper is organized as follows. In Section|2]we describe 
our model. In Section[3]we present our algorithm for finding 
best fitting models. We present our results for four dSphs in 
Section We discuss the results in Section [5] and give our 
conclusions in Section|S] 

2. MODEL 

2.1, Overview and Assumptions 

We assume that initially gas was isothermal and in hydro- 
static equilibrium in the static potential of the DM halo. Stars 
are assumed to have formed from gas instantaneously with a 
certain efficiency £ out to a radius where the gas density drops 
below some critical value. The velocity dispersion of newly 
born stars is assumed to be equal to the sound speed in the 
gas. After the leftover gas is expelled instantaneously by the 
combined action of ionizing radiation, stellar winds, and su- 
pernova explosions, the stellar cluster is relaxed over a few 
crossing times through violent relaxation and/or phase mix- 
ing. We simulate this process numerically using an A^-body 
code. 

As we will show in this paper, in our model the stellar clus- 
ter does not expand much after the star burst, resulting in the 
final central velocity dispersion being only ~ 20% smaller 
than the initial value. The central line-of-sight velocity disper- 
sion f or dSphs is in the narrow interval 8 .5 ± 2 km s" 1 ( Mateo 
1998). It is remarkable that this parameter is almost constant 
for all Galactic dSphs. It is also comparable to the sound 
speed in the warm neutral mediu m, 6—7 .5 km s" 1 , in the ISM 
of the Milky Way ( Wolfire et al. 1995; we used their "stan- 
dard" model and assumed that the mean molecular weight is 
14/1 lm p , where m p is the mass of the proton), and ISM ve- 
locity dispersion of ~ 9 ± 2 km s" 1 in dwarf irregular galaxies 
(Mate oll998l) . The vertical velocit y dispersion of young stars 
in the Galactic disk, 6 ± 3 km s" 1 ( Haywood, Robin, & Creze 
is also in the same range, giving support to our model 
assumption that the velocity dispersion of newly born stars is 
comparable to the sound speed in the star-forming gas. 

To justify our assumption of the gas being isothermal, we 
should mention that in all of our models the relaxed stel- 
lar clusters are almost isothermal (in terms of total stel- 
lar velocity dispersion). Galactic dSphs are known to have 
roughly isothermal line-o f-sight velocity dispersion p rofiles 
in their central regions ( Mate3 H997l iKlevna et al] l200l 
iWilkinson et alJl2004l) . We can also reproduce the observed 
decline in the line-of-sig ht velocity dispersion in the outskirts 
of some dSphs (Kleynaet al. 2004), which in our model is 
caused by a fraction of stars which move beyond the initial 
stellar cluster radius during the dynamical relaxation phase, 
and as a result have strong radial anisotropy. 

Because the DM halo potential in our model is static, nei- 
ther adiabatic contraction of the central part of DM halo 
due to presence of baryons nor adiabatic expansion after 
the removal of the remaining gas in a post-starburst sys- 
tem are modeled. Static potential models presented in this 
paper took around 500 CPU-days to compute; live DM 
halo simulations would take orders of magnitude longer to 



run, which is obviously not feasible. From both observa- 
tional and theoretical points of view, the situation with in- 
ner density profiles of DM halos is less than clear. Obser- 
vations produce less steep inner density slopes for galactic 
DM halos (7 = -1 . . .0) than the predictions of cosmologi- 
cal CDM models (7 = - 1.5 ... - 1). Though a few me cha- 
nisms [such as warm DM (Bode, Ostriker, & Turok 2001), re- 
coiling black holes (Bo vlan-Kolchin. Ma. & Ouataerl 2004). 
and dynamic heating due to presence of dense gas clumps 
( El -Zant, Shlosman, & Hoffman 2001)] were proposed to 
bring theory in agreement with observations, none of these 
explanations has received yet wide acceptance. Given the un- 
certainty in the inner density slope for galactic DM halos, we 
believe that by considering two different types of static DM 
halos, with flat core (7 = 0) and with 7 = -1 cusp, we allevi- 
ated some inadequacies resulting from the use of static DM 
potential. 

2.2. Initial Gas Distribution 

Hydrostatic equilibrium of self-gravitating gas inside the 
static potential of a spherical DM halo can be described us- 
ing the following expression: 

1 dPg _ ggtot (1) 
p g dr dr 

Here P g and p g are the pressure and the density of gas, and 
<E> tot is the total gravitational potential with contributions from 
both gas and DM. The radial gradient of gravitational poten- 
tial in equation Q can be calculated as d$ tot /dr = G[M(r) + 
m g (r)]/r 2 , where M(r) and m g (r) are the enclosed masses of 
DM and gas, respectively, and G is the gravitational constant. 
To close equation Q, one has to write Poisson equation for 
the gas component: 

Here ip g (r) is the gas contribution to the total potential. 

W e consid e r two types of DM halo density profiles: flat- 
core iBurkerll i ll 9951) halos and Navarro-Frenk- White (NFW) 
halos which have a central d ensity cusp with 7 = -1 
iNavarro. Frenk. & White! Il997l) . The density profiles for 
these halos are as follows: 

M- axm*m ,Bmkml m 

M' TjJhw (NFW) - (4) 

Here r s and po are the scaling radius and density for DM halos. 
(In Burkert halos, po corresponds to the central DM density; 
in NFW halos, which have divergent central density, p = po at 
~0.5r s .) 

For the following analysis, it is convenient to switch to 
new dimensionless variables: x = r/r s for radial distance, 
p' = p/po for density, M' = M/(47rr^po) for mass, v' = 
v/(4TrGijp Q y/ 2 for velocity, and t ' = t (4ttGpq) 1/2 for time. In 
these new variables, the enclosed DM mass for the two types 
of halos can be written as follows: 

r (1 /2)[ln(l +x) + (l/2)ln(l +x 2 )- 
M'(x)=< -arctanx] (Burkert), (5) 

Un(l+x)-.x/(l+x) (NFW). 
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We assume that gas is isothermal: P g = p g ci (here c g is the 
sound speed in gas, which is assumed to be constant). Equa- 
tions Q and (0 then can be rewritten as 

d A = -A-^ [M '( X ) +m > g i (6 ) 

ax x L s 

and 



Equations (|5j{7) form a system of two ordinary differential 
equations for two variables, p' g (x) and m' g (x), which has to be 
solved numerically. The boundary conditions are p' g (0) = p' g0 
(a free parameter) and m' g (0) = 0. 

2.3. Star Formation and Dynamical Relaxation of Stars 

We adopt a simple star formation model. We assume that 
a certain mass fraction £ of gas is turned instantaneously into 
stars everywhere inside the halo where the dimensionless gas 
density is above a certain threshold value A (in units of the 
central gas density p' ?0 ). In other words, star formation effi- 
ciency is assumed to be £ (a constant) for p' g > Xp' g , and zero 
otherwise. The parameter A was introduced to avoid the star 
formation occuring at unphysically low densities. The newly 
born stars are assumed to have an isotropic distribution of ve- 
locity vectors, with the one-dimensional velocity dispersion 
equal to the sound speed in gas c'. 

After the starburst, all the remaining gas is assumed to be 
removed from the system instantaneously by the combined 
action of ionizing radiation, stellar winds, and supernovae ex- 
plosions. In reality, a fraction of this gas can be later reac- 
creted by the halo leading to subsequent episodes of star for- 
mation. For simplicity, here we consider all stars in a dwarf 
spheroidal galaxy to be formed in a single starburst. 

It is easy to see that models with £ = 1 and A = will re- 
main in equilibrium after all gas is turned into stars. Indeed, 
our hydrostatic gas equilibrium equation Q is identical to the 
isotropic stellar Jeans equation if = p g and = c g . (Here 
and er* are density and one-dimensional velocity dispersion 
of stars.) 

For the cases with £ < 1 and/or A > 0, the initial configura- 
tion of stars after the gas removal is not an equilibrium one: 
stars are dynamically hot, and will expand before reaching 
an equilibrium state. This highly non-linear process involves 
phase mixing and/or violent relaxation. Consequently, we had 
to resort to numerical A^-body simulations to derive the final 
equilibrium configuration of our stellar clusters. To run the 
simulations, we used the pa rallel version of the p opular mul- 
tistepping tree code Gadget (Sprin g^! et alJ20 01). 

2.4. Numerical Scheme 

Overall, there are four free parameters in our model which 
can influence properties of relaxed stellar clusters: the sound 
speed in gas in units of the halo scaling velocity c' g , the central 
density of gas in units of the halo scaling density p' g , the star 
formation efficiency £, and the minimum gas density A for 
star formation to take place in units of the central gas density. 

The parameter c' can be considered as a measure of how 
hot gas is in relation to the virial temperature of the halo. 
The nuance is that in our model the gas temperature is mea- 
sured in terms of "scaling temperature" which can be defined 



as T s = pVf /(2k), whereas the usual definition for the virial 
temperature of a halo is r v i,- = fiV^/flk). (Here p is the mean 
molecular weight of the gas, k is the Boltzmann constant, V s 
and Vyir is the halo circular speed at the scaling radius r s and 
the virial radius r vlr , respectively.) The connection between 
T s and T vn depends on halo density profile (Burkert or NFW) 
and concentration c = r w \ r /r s . For example, for both Burkert 
and NFW halos with c = 3.5 ... 10, gas in our models is at the 
virial temperature when c' g = 0.46 . . . 0.39. 

We consider three values of the parameter A: (no star 
formation cutoff), 0.1, and 0.3. We cover a wide range of 
values for the remaining three parameters: c' g = 0.01 . . .0.48, 

/^ )0 = 0.1...10 5 ,and£ = lCT 3 ...l. 

To generate an initial particle distribution for the A^-body 
simulations, we use the enclosed gas mass function miix) 
which is derived numerically by solving equations (|5jQ. 
First, we find the star formation cutoff radius x\ such as 
p' g {x\) = A/5g . (For the case of A = 0, x\ is set to be 
equal to the largest considered radius x max = 10 3 ). Next, the 
integral probability function P(x) is constructed as P(x) = 
m' g (x)/m' g (x\). For each A^-body particle, the radial dis- 
tance x is obtained by generating a uniform random num- 
ber y = [0...1] and then solving numerically the equation 
P(x) = y. Finally, to obtain three components of the radius 
and velocity vectors we assume that they are randomly ori- 
ented in space, which is true for a spherically symmetric sys- 
tem with isotropic distribution of velocity vectors. We as- 
sumed Maxwellian distribution of velocities, with the one- 
dimensional velocity dispersion equal to the sound speed in 
the gas c'. 

All our models were simulated with N = IQ 4 equal mass 
stellar particles. This number is large enough to allow for 
accurate surface brightness profile calculations over a large 
range of radial distances and to avoid significant spurious dy- 
namical evolution caused by close encounters between parti- 
cles. At the same time, N is small enough to make it feasible 
to simulate hundreds of models. The masses of individual 
particles are equal to m*/N, where = £m g {x\) is the total 
stellar mass. 

The value of the gravitational softening length e was cal- 
culated as Q.llx h N~ lli JHavashi et alJ l2003). where the ini- 
tial stellar half-mass radius Xf, was derived by solving the 
equation P(x/,) = 1 /2. The values of the Gadget parame- 
te rs which control accuracy of in tegration were the same as 
in Mashchenko & Sills (2005a b), where a total energy con- 
servation of better than 2% was achieved when simulating a 
globular cluster for 6000 crossing times at the half-mass ra- 
dius. In particular, the individual time steps in the simula- 
tions were equal to (2r]e/a) x l 2 , with the accuracy parameter 
77 = 0.0025. (Here a is the acceleration of a particle.) The 
parameter controlling the force calculation accuracy (EnTol- 
ForceAcc in Gadget notation) was set to 0.01. 

In the A^-body simulations, the DM halo is represented by a 
static gravitational potential, which is identical to the one used 
for solving equations ||5J0. The important point to make is 
that our DM halos are not truncated at the virial radius — 
neither in solving equations nor m A^-body simulations. 
We did the consistency check (that the observed size of the 
galaxy is smaller than the inferred virial and tidal radii) after 
fitting the models to observations. 

We ran our A^-body models for the interval of time t' = 370 
(in dimensionless units), with 200 intermediate snapshots per 
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model. This time interval is much larger than the crossing 
time f^ ross = 1 jc' at the scaling radius r s , which is ^ 100 for 
all our models. 

To cover well the 4-dimensional free parameter space, we 
had to run a significant number of models (around 700). To 
facilitate the process, we designed a pipeline for automatic 
generation of observable properties of the models for given 
DM halo profile (Burkert or NFW) and given values of the 
four initial free parameters c' g , p' , £, and A. The pipeline 
consists of the following steps. 

1 . System of differential equations (|5j0 is solved numeri- 
cally, producing isothermal hydrostatic gas distribution 
in the potential of the DM halo. 

2. Using the procedure described above in this section, the 
initial distribution of stellar particles is generated. 

3. An Af-body code is run to simulate the stellar cluster 
relaxation in the static DM halo potential. 

4. For last Ns = 100 snapshots (corresponding to t' = 
185... 370), both time-averaged values and standard 
deviations from mean are calculated for surface bright- 
ness E in different radial bins. To calculate E, we 
us e the projection method d escribed in Appendix B 
of Mashchenko & Sills! (I2005al) . The noisier innermost 
and outermost parts of the surface brightness profile are 
cut off at the radii where the standard deviation from 
mean becomes larger than A = 0.2 dex. In this fash- 
ion, we ensure that the inner part of the E profile has 

1 /2 

converged to better than A/N s = 0.02 dex (because 
of very short crossing time near the halo center result- 
ing in virtually uncorrected surface brightness values 
in different snapshots), and that the outer part of the 
profile (where a strong correlation is expected) has con- 
verged to 0.02 . . .0.2 dex. We repeat this procedure for 
a smaller number of snapshots Ns = 50 (corresponding 
to t' ~ 280. . .370), and choose the profile which has 
converged over a larger range of radial distances. 

5. We calculate a line-of-sight velocity dispersion <jq 
[usin g a procedure delinea ted in Appendix B of 
Mashchenko & Sills (2005a)] averaged over the core 
region of the galaxy (within the half -brightness radius) 
and time-averaged over the last Ns snapshots. 

As a test, we compared the radial stellar density profiles 
in our evolved models with £ = 1 and A = with the corre- 
sponding initial gas density profiles (they should be identical 
— see § 12.31 . In most of the models, the profiles are found 
to be identical within the measurements errors. In a few high 
density models, with p' = 10 3 and 10 5 , we see spurious flat- 
tening of the inner stellar density profile due to interaction 
between stellar particles. This effect is found to be negligible 
at the radii larger than the softening length e. 

2.5. General Properties of the Models 

Surface brightness profiles in our relaxed models exhibit a 
variety of shapes — from simple cases of monotonously in- 
creasing slope, to profiles with one or more inflection points. 
The important fact is that in most cases we observe a flat 
core plus extended halo (with different slopes) in the surface 
brightness profiles for both Burkert and NFW static halos — 
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FIG. 1. — Dependence of the estimated radial expansion factor xlxm\ on 
initial central gas density p' Q (in the limit £ — > 0). Thick and thin lines 

correspond to models with = 0.16 and 0.01, respectively. Solid (dashed) 
lines correspond to NFW (Burkert) halos. In all the models A = 0. 



somewhat like profiles of dSphs. An obvious central stellar 
cusp is seen only in the models with strongly self-gravitating 
gas: p' g > 10 3 (for both Burkert and NFW halos). 

Our relaxed stellar models are approximately isothermal 
(in terms of total three-dimensional velocity dispersion). For 
models with non-zero star formation lower density cutoff, A, 
the line-of-sight velocity dispersion profiles are often sharply 
declining at large radii. For such models, the outer stellar halo 
is strongly radially anisotropic. 

On average, our relaxing stellar clusters do not expand 
much, and as a consequence do not cool down much adiabati- 
cally. For all our models, the ratio of initial (non-equilibrium) 
core line-of-sight velocity dispersion cr^o to its final, relaxed 
value do varies between -0.01 and 0.43 dex, with a small av- 
erage value of 0.10 ±0.07 dex. To explain this model prop- 
erty, let us consider the extreme case of £ — > resulting in the 
largest possible expansion factor. Initially, stars are born on 
equilibrium orbits (due to the fact that our hydrostatic equilib- 
rium equation is ident ical to the isotropic stellar Jeans equa- 
tion - see Section l2~3l . An instantaneous removal of virtually 
all baryons after the starburst does not change individual or- 
bital angular momenta of the stars, resulting in the following 
adiabatic law: 



a 



(8) 



Here x is the baryonic half-mass radius and a' is stellar ve- 
locity dispersion. Suffix "ini" stands for the initial, non- 
equilibrium values. The virial equation for the relaxed cluster 
in case of £ — > is 



M'(x) 



x M'(x w ) + m'/2' 



(9) 



where the gravitational radius is assumed to be equal to the 
half-mass radius. We estimate the largest expansion factor 
x/xi n i by combining equations (jSJl and (|9ji: 
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FIG. 2. — Radial profiles for Burkert models with c' g =0.16, p' ? = 10, and § = 0.1. (a) Stellar density profiles. Thick solid line shows the initial non-equilibrium 
profile. Profiles for relaxed models with A = 0, 0.1, and 0.3 are shown as thin lines (solid, long-dashed, and short-dashed, respectively). The dotted line shows 
the DM density profile, (b) Stellar one-dimensional velocity dispersion cr' (left ordinate) and anisotropy j3 = 1 — (c'/cr') 2 (right ordinate) as a function of radial 
distance x for model with A = 0. 1. Thick solid, short-dashed, and long-dashed lines correspond to total dispersion cr', radial dispersion cr', and tangential 
dispersion cr,' for the relaxed model. Horizontal dotted line shows the initial non-equilibrium value of velocity dispersion. Vertical dash-dotted line marks the 
cutoff radius x\ . Thin solid line shows radial /3 profile for the relaxed model. 



x ^ M'(x mi ) + m' g /2 
x ini ~~ M'(x) 



(10) 



TABLE 1 

Observational data on dwarf spheroidals 



Assuming that most of the gas is concentrated within the scal- 
ing radius of the DM halo (which is the case for almost all 
of our models), M'(x) can be approximated as M'(x) oc x a , 
where a = 2 for NFW and 3 for Burkert halos. In this case, 
the expansion factor can be written as 



1 + - 



2M'(x ioi ) 



l/(a+l) 



(ID 



It happens, that for isothermal gas in hydrostatic equilibrium 
inside a static DM halo, the estimated expansion factor x/x- m i 
is never larger than ~ 1 .7 for the whole range of considered 
c' g and p' , and for both NFW and Burkert halo density pro- 
files (see Figure^. For less extreme cases (with larger £) we 
expect the expansion factor to be even closer to unity. 

To illustrate the impact of a lower density cutoff for star 
formation A on our models, we show in Figure |2 a few ra- 
dial profiles for Burkert models with c' = 0.16, p' = 10, and 
£ = 0.1. From the panel (a) of this figure one can see that non- 
zero values of A both change global properties of the relaxed 
stellar clusters (such as central density and size) and modify 
the shape of the density profile. In panel (b) of Figure|2]the ra- 
dial velocity dispersion profiles are shown for the model with 
A = 0. 1 . (This model is close to some of the best fitting mod- 
els for Galactic dSphs — see Section^) As you can see, stars 
in the relaxed cluster show strong radial anisotropy (with the 
anisotropy parameter j3 — * 1) outside of the initial radius of 
the cluster x\ (the vertical dash-dotted line). The total veloc- 
ity dispersion (thick solid line) does not change much with 
radius. The transition from an almost isotropic to a strongly 
radially anisotropic distribution of velocities is very sharp in 
this model. It is interesting to note that in the outskirts of 



Name 
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So 


°"obs,0 


T 






kpc 


mag arcsec~ 2 


km s 






Draco 


82 


25.3 


9.5 ±1.6 


1.32 


13.2 a , 147.0 b 


Sculptor 


79 


23.7 


6.6 ±0.7 


1.23 


296.1 


Carina 


101 


25.5 


6.8±1.6 


0.85 


28.7 


Fornax 


138 


23.4 


10.5 ± 1.5 


0.94 


304.9 



NOTE. — Here XKing ' s me X 2 value for the best fitting theoretical King 

model. 

a For the dataset Draco 1 [from Odenkirchen et al. COO 11) 1. 
b For the dataset Draco2 [from Wilkinson et al. 1 2004|)]. 



the cluster the radial velocity dispersion (short-dashed thick 
line) is comparable to the initial velocity dispersion (horizon- 
tal dotted line). 

3. COMPARISON WITH OBSERVATIONS 

3.1. General Remarks 

We use two types of observational data on Galactic 
dwarf spheroidal galaxies in our analysis: star number den- 
sity profiles (with corresponding one-sigma uncertainties) 
and co r e line- of-sight velocity dispersion (T b s .o [taken from 
Mateo (1998), and listed in our Table The star num- 
ber density profiles for four dwarf spheroidals analyzed in 
this paper (Dra co, Sculptor, Carina, and Fornax) are taken 
from Odenkirc hen et alJ feOOll) . IWilkinson elaH ( I2004 . and 
Walcher et al. (2003). For our dwarf galaxies, we adopt dis- 
tance Z) and central surface brightness So values from Mateo 
(1998), and baryon mass-to-light ratio T from iMateo et al 
( 1998, their Table 6; we averaged T values for Salpeter and 
"Composite" initial mass functions). These three parameters 
(D, So, and T; listed in Tabled are use d to convert the ob- 
served star number density profiles into stellar surface density 
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profiles E(r) (in M pc~ 2 units). To perform the conversion, 
we estimated the central star number density in the observed 
profiles by averaging values for a few central radial bins, well 
within the galactic core. 

Due to the time-consuming nature of the simulations, high 
dimensionality of our free parameter space, and some degree 
of noisiness in the results of our simulations, we do not at- 
tempt to find a model which is the best match to observa- 
tions. Instead, for each galaxy we try to find at least one (typ- 
ically more) models which would simultaneously satisfy the 
following conditions: (1) Surface brightness profile shou ld be 
very similar to the observed one (quantified in § 13 .21 . (2) 
Core (within the half-brightness radius) line-of-sight veloc- 
ity dispersion <tq should be equal to the observed value <7 bs,o 
(Table within the measurement errors. (3) The values of 
the DM halo scaling mass M s = Airrlpo and scaling radius r s , 
which are inferred after fitting the model surface brightness 
profile to the observed one, should correspond to realistic ha- 
los pred icted by cosmological ACDM simulations (quantified 
in §13). 

3.2. \ 2 Fitting of Surface Brightness Profiles 



To quantify the first of the selection criteria in § 13.11 we 
measure weighted \ 2 between the model and observed sur- 
face brightness profiles (in linear space). We minimize \ 2 nu- 
merically as a function of two parameters — a multiplier 5s 
to transform the model column density to the observed one, 
and a multiplier S r to rescale the linear size of the model. 
Other rescaling coefficients can be derived from the above 
two parameters: S p = Ss/S r for density, S„, = S^S 2 . for mass, 
S v = (SsS,-) 1 / 2 for velocity, Sj = SsS,- for initial gas tempera- 
ture, etc. For example, we can derive the scaling mass M s and 
scaling radius r s for the DM halo. 

The minimum \ 2 values are much larger than unity for all 
four galaxies (this is also true for best-fitting King models; 
see Table 0. The reasons for that are that (1) the observed 
surface brightness profiles show many fine details which can- 
not be reproduced by any simple model, and/or (2) the quoted 
measurement errors are underestimated. As a result, our strat- 
egy is to look for models producing smaller \ 2 values than the 
best fitting King models: x' 2 = X 2 lx\ 



3.3. Plausibility of DM halo s 

To quanti fy the third cri terion in 33.11 we used the formula 
of Sheth & Tormen ( 1999) to calculate the comoving number 
density of DM halos as a function of virial mass M v ; r and red- 
shift Z- 



dn 



d\nM v 



1 + - 



PmV 

s 



dS 



dM v 



exp 



y ) • d2) 



Here v = (a/S)'/ 2 <5(z), the present day mass density of the 
universe is p m = 3f2 m // 2 /(87rG), and the values of numerical 
coefficients are A = 0.322, a = 0.707, q = 0.3; fl m and H are 
mass density of the universe in critical density units and Hub- 
ble constant. We calculate the variance S of the primordial 
density field on mass scale M v ; r , extrapolated l inearly to z = 0, 
from fi tting formulae given in Appendix A of van den Bosch 
(2002), which are accurate to better than 0.5% over the mass 
range 10 6 < M vrr < 10 16 M Q . The shape parameter T, re- 
quired for the calculation of S, is estimated as 



(13) 



r = n m hexp[-n b (i + V2h/n,„)] 



JSugivamall995l) . where h = #/(100 km s" 1 Mpc -1 ) and Q. b 
is the baryonic mass density of the universe in critical den- 
sity units. The critical overdensity threshold for spherical col- 
lapse, linearly extrapolated to z = 0, is approximated by 

S(z) = 0.15(12 7 r) 2 / 3 ^ 1 0055 /£>(z) (14) 

JNavarro et all 19971) . which is valid for flat ACDM cosmolo- 
gies. The fitting formula f or the linear growth factor D(z) 
was taken from Appendix of iNavarro et alJ Through- 
out this paper, we adopt the following values of cosmological 
constants: fl m = 0.27, Q b = 0. 044, H = 71 km s" 1 Mpc" 1 , and 
cr 8 = 0.84 ( Soer gel et all2003l) . We assume a flat ACDM cos- 
mology: £1a = 1 — fl m . 

Fitting model surface brightness profiles to the observed 
ones gives us the values of t he s caling halo mass M s and 
scaling halo radius r s (see § 13. 21 . To use equation d!2l >. 
we need to know M v ; r and z. To make a transformation 
from (M s ,r s ) to (M v ; r ,z), we use the fact that in cosmologi- 
cal ACDM simulations the halo concentration c = r v ; r / r s dis- 
tribution is approximatel y lognormal (for fix ed M v ; r and z), 
with dispersion 0.14 dex ( Bu llock et all 2001). For low mass 
halos (with M v ; r = 10 8 ... 10 11 Mq), hal o co ncentrati on de- 
pends on Myj r and z in the following way ( Bulloc k"et all2001t 
ISternberg etalJ2002l) : 



27 

c=— 10' 
1+z 



0.14^ 



M v 



(io-a/. 



-0.08 



(15) 



Here v c is the number of standard deviations from the 
mean concentration. Multiplying dn/ dh\M wu (eq. 1121 ) 
by the probability to find a halo with v' c = v c ...v c + dv c , 
(2ir)~ x l 2 exp(-Vc /2)dv c , and dividing the result by dv c gives 
us the comoving halo number density per unit lnM v j r and per 
standard deviation in concentration: 



dn 



d\nM w ivdv c 



A 
2^ 



1 + - 



PmV 
S 



dS 



x exp 



t/Mvir 

v 2 + v} 



(16) 



We use equation d!6t to find the most probable combination 
of M v ; r and z for a halo with given M s and r s . Our algorithm is 
as follows. We vary v c in the interval -4 ... 4 with a small step 
(0.01). For each value of v c , the following non-linear equation 
is solved numerically to find z: 



Here 



Mvir(z)=M s M'{c[Mvir(z),z]}. 



Mvir(z) = [(27r s lO ol4l/ O 3 (lO 9 M ) a24 x 



(17) 



A r H 2 n„ 



-, 1/(1+0.24) 



2GQ. 



(18) 



was derived from the definition of critical overdensity for col- 
lapsed halos 



2GM v i 1 U„. z 
rl.H 2 n m {l+z) 3 



(19) 



and equation J15i for halo concentration c; M'(x) is given by 
equation (|5j- The mass density in flat universe as a function 
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of redshift is O mz = {1 + Q A /[O m (l + z) 3 ]} We use the fit- 
ting fo rmula A e ~ 187r 2 + 82<i-39<i 2 from lBrvan & Normanl 
( 1998), which is accurate to 1% in the range Q mz = 0.1 - 1. 
(Here d = fi OT)Z — 1.) 

After solving equation (I17> (which gives us z for given M s , 
r s , and v c ), we use equation Jl 81 to find the corresponding M v ; r 
value. Finally, from equation (I16> we estimate the comoving 
halo number density F(v c ). Repeating the above procedure 
for the whole interval of v c , we find the v c value (and corre- 
sponding values for z and M vlr ) which maximizes the likeli- 
hood function F(v c ). 

We make the following very rough estimate of the lower 
cutoff F m i n value (such as that the halos with F <C Fmin 
would be very unlikely progenitors of dSphs, whereas the ha- 
los with F > Fmin would be likely progenitors). We divide 
the esti mated mass of the Local Group of ~ 2.3 x 10 12 Mq 
ivan den Berg hll999) by the present day mass density of the 
universe p„, to derive the comoving volume of the part of 
the early universe which became the Local Group, Vlg ~ 
61 Mpc . (This volume corresponds to a sphere with comov- 
ing radius of 2.4 Mpc.) There are A^jsph ~ 20 dSph galax- 
ies in the Local Group. From the analysis presented in this 
paper, the DM halo masses for these objects span ~ 2 dex, 
so AlnMvk ~ 4.6. When the very first of these halos was 
virialized, the corresponding comoving number density was 
F ~ l/(VL G AlnM vir ). Much later on, when all JVds P h ha- 
los had formed, the comoving number density became F ~ 
NdSph/(VLGAlnMvi r ) (which would only be accurate if the ha- 
los did not accrete any mass after their formation). To come 
up with a single number, we take an average (in log) of the two 

above extremes: (F) ~ A^/(7ujAlnMvir) ~ 0.016 Mpc -3 . 

We adopt F m ; n = 0.01 Mpc" 3 for the rest of this paper. 

3.4. Strategy of Searching for Best Fitting Models 

In this project we faced the formidable task of sampling 4- 
dimensional (c g , p' g0 , £, and A) initial parameter space in our 
numerical simulations (with an extra factor of 2 due to the two 
different types of DM halo profiles considered — NFW and 
Burkert). As we mentioned before, the requirement for gas to 
be colder than the virial temperature results in < c' g < 0.5. 
Other obvious constraints are < £ ^ 1 and ^ A < 1 (by 
definition). We had no other reasonably motivated constraints 
on the range of our free parameters. 

Initially, we sampled a large range of the free 
parameter space by simulating a grid of models 
with c' g = (0.01,0.04,0.16), p' g0 = (0.1,10,10 3 ,10 5 ), 
£ = (0.001,0.01,0.1,1), and A = (0,0.1), which resulted 
in 192 models (both NFW and Burkert). We quickly realized 
that models with c' g < 0.04 or p ' q > 10 can be safely 

ruled out as they all produced unacceptably large x 2 values 
(in other words — their surface brightness profiles were 
nothing like profiles of dSphs) and/or very low F values 
being much smaller than F m m (meaning that the inferred DM 
halo parameters were improbable in the adopted ACDM 
cosmology) for our sample of galaxies. 

Other important observations were: (1) the range for c' g 
should be expanded up to ~ 0.5 to include warmer models 
(which we could not do for models with A = 0, as in most of 
such models the initial stellar half-mass radius becomes diver- 
gent for c' g > 0.2); (2) in the range 0.04. . .0.5, the parameter 
c' should be better sampled (by a factor of two) ; (3) in the 
range 0.1... 10 3 , p' g should also be better sampled (by a fac- 



tor of two). 

As a result, we ran 294 additional models, which also in- 
cluded a set of models with A = 0.3 (for such models we 
simulated only the cases with c' g ^ 0.04 and p' gfi < 10 3 ). 
Altogether, we had 486 models on a regularly spaced 4- 
dimensional grid. 

The problem we faced at that point was that no realistic 
number of models on a regularly sampled 4-dimensional grid 
would be large enough to simult aneo usly satisfy all our three 
selection criteria (outlined in § 13 . II . More specifically, de- 
spite the fact that for most of our galaxies a relatively large 
volume in the 4-dimensional parameter sp ace would satisfy 
both the first and the third criteria from § 13.11 few (usually 
none) of the models from this volume would produce the core 
line-of-sight velocity dispersion <jq which would be consistent 
with observations within the measurement errors (our second 
criterion). 

Fortunately, we found a solution for the above problem. We 
noticed that for the models which produce reasonably good \ 2 
fits (within a factor of two from the model with the lowest % 2 
value), there is a good multivariate correlation (in logarithmic 
space) between oq and the free parameters c' g , p' g0 , and £ for 
NFW models, and p' g and £ for Burkert models (the parame- 
ter A is fixed). 

Our strategy was then as follows. For each galaxy, each 
value of A, and each type of DM profile (NFW or Burkert), 
we used a multivariate \ 2 fitting to find coefficients for the 
ao(c g ,p' ,£) correlation. (Only the models producing good 

surface brightness \ 2 fits were used.) We used this correlation 
to randomly sample a plane <7 bs.o = co(Cg,PgO'0> with typi- 
cally 10 additional models for each case. (Here er b s ,o is the 
observed core line-of-sight velocity dispersion, listed in Ta- 
ble ^) The sampling was done within the three-dimensional 
rectangular volume corresponding to models which simulta- 
neously satisfied our first and third criteria from § 13.11 As 
our method is not precise, not all additional models satisfied 
all the three criteria. Figure[3]illustrates our algorithm for the 
case of the Draco 1 dataset (see Section 0, NFW halos, and 
the A = 0.1. 

Overall, we simulated around 700 models. It took approx- 
imately 500 CPU-days to simulate all the models on McKen- 
zie cluster at CITA composed of 256 dual Intel Xeon 2.4 
GHz nodes networked with commodity gigabit ethernet. Each 
model was run on 4 CPUs, with many models running in par- 
allel. 

4. RESULTS FOR FOUR DWARF SPHEROIDAL GALAXIES 

4. 1 . General Remarks 

In Table |2] we show statistics on " good " models (which 
meet all three selection criteria from § 13 . 1 1 for five datasets 
considered in this paper — for the Draco (two datasets), 
Sculptor, Carina, and Fornax dSphs. Because of the small 
number of "good" models N g , the not completely random na- 
ture of our best fitting models search algorithm, and some de- 
gree of noisiness in the results of our simulations, in this table 
we present three complimentary statistical measures for each 
parameter and for each dataset: the value for the "best" model 
(producing the smallest value of \ 2 m surface brightness pro- 
file fitting), the mean value with the associated standard de- 
viation, and the total range for all N g models. Most of the 
parameters are presented as their logarithms as their distribu- 
tion is much less skewed in logarithmic space. For the same 
reason, instead of the redshift, z, of halo formation, we present 
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FIG. 3. — Illustration of our best fitting model search algorithm for the Dracol dataset (see Section^), NFW halos and A = 0.1. Panels (a-c) show the projection 



of the models onto three orthogonal planes in the three-dimensional space of free model parameters 



g< r g ,o 



and f , Crosses mark models which have unphysical 



DM halos (F < F m [ n ). Small open circles correspond to physical (F > F m i n ) models which have a bad x 2 fit between the modeled and observed surface brightness 
profiles (x 2 > 2Xmin' where Xmin ' s tne lowest \ 2 value for all ~ 700 models). Solid circles show the location of models which both are physical (F > F m i R ) and 
fit well the observed surface brightness profile (x 2 < 2x 2 nin )- The large circle marks a physical model with the smallest \ 2 value (for NFW halo and A = 0.1). 
Rectangular boxes delineate the area where good models are being looked for. Panel (d) shows the multivariate correlation between the three free parameters (c 1 ,, 
p' g , and §) and the core line-of sight velocity dispersion uq for models with x 2 < 2Xmj„. which is used for searching for good models. Here crosses correspond 



to unphysical DM halos (with F < F ra i„), whereas the solid circles mark physical halos (F > F m j„). Two horizontal dotted lines show the one-sigma errorbars for 
the observed <T o b s .0 in Draco. 



statistics on the cosmic scale factor a = 1/(1+ z). 

We use the foll owing definition of tidal radius rta 
JHavashi etalJ20 03): 



M( m ) 



'tid 



R dM H 



M H (R) dR 



M H (R) 
Ri ' 



(20) 



Here M(r) and M H (R) are enclosed mass functions for the 
satellite and host galaxies. We assumed that Milky Way DM 
halo has NFW density profile with M vil = 1.5 x 10 12 M Q 
(corresponding to R v ; r = 298 kpc). From equation (14) of 
ISternberg et all J2002I) . the concentration of the Milky Way 
halo is c = 13.2. The tidal radius was calculated for the current 
distances D (listed in Tabled of the dSphs from the Galactic 
center. 

To calculate the temperature T of the gas, we assumed that 
the mean molecular weight is \l = 1 .218m p , which is appropri- 
ate for neutral gas of primordial composition. (Here m p is the 
mass of the proton.) If gas is ionized, the temperature would 
be a factor of 2 smaller. 

In the next four sections we also present surface density 



E(r) and line-of-sight velocity dispersion a(r) profi les for 
our best fitting models for each dataset (Figures |4I8> . Only 
those parts of the profiles which have a standard deviation 
(estimated from the last Ns snapshots for a given radial bin) 
smaller than 0.2 dex for E(r) and 1 km s" 1 for <j{r) are shown. 
For comparison, we s how the surf ace density profiles for the 
best fitting theoretical King ( 1966) models (dashed lines). 

4.2. Draco 

Draco is probably the best candidate for having an extended 
DM ha lo among Galactic dSph satellites. Odenk irchen et alJ 
(2001) used deep multicolor photometry covering 27 square 
degrees around the center of this galaxy to show that Draco 
has a very regular structure (with no signs of tidal interaction 
with the Milky Way gravitational field) out to the radius where 
its surface brightness falls to 0.003 of the central value. Their 
results suggested that there is much more DM in the outskirts 
of Draco than in its center. 

More recently, Wilkinson et al. (2004) presented a new star 
number count profile for Draco. The new profile appears to 
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FIG. 4. — Best fitting models for the Dracol dataset. (a) Surface density profiles. The observed profile is shown as solid circles with errorbars. The solid 
line corresponds to the best fitting model (with the lowest x' 2 )- The dashed line shows the best fitting theoretical King model profile, (b) Line-of-sight velocity 
dispersion profiles for all 1 1 "good" models from Tablel2l The thick line corresponds to the model with the lowest x' 2 ■ 
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FIG. 5. — Best fitting models for the Draco2 dataset. (a) Surface density profiles. The observed profile is shown as solid circles with errorbars. The solid 
line corresponds to the best fitting model (with the lowest x' 2 )- The dashed line shows the best fitting theoretical King model profile, (b) Line-of-sight velocity 
dispersion profiles for all 13 "good" models from Tablel2l The thick line corresponds to the model with the lowest x' 2 ■ 



be even more shall ow in the outskirts of the galaxy than the 
profile of Odenkirchen et al. (2001 ), though the tw o profiles 
are marginally consistent. Wilkinson et al. (2004) also pub- 
lished a line-of-sight velocity dispersion profile for Draco, 
which is roughly isothermal, with the exception of the last 
radial bin where the stars are marginally colder than the rest 
of the galaxy. 

We use star count profiles from both lOdenkirchen et al] 
d200lL their sample S2, which has a lower estimated fore- 
ground/background contamination level than their sample SI) 
and Wilkinson et al. (2004) in our analysis, which we desig- 
nate Dracol and Draco2, respectively. We found 11 and 13, 
respectively, "good" models for these two datasets (see Ta- 
ble |2jl. In both cases, our models produce much better \ 2 



surface brightness fits than best fitting King models (x = 
xVxking ~ 0.6). As can be seen in Figures^ and|5^, the 
only statistically significant difference between our best fitting 
models and King models is in the outer parts of the galaxy. In 
the Draco outskirts, the King models exhibit a sharp cutoff 
due to the external tidal field, whereas our models have a rela- 
tively shallow power-law profile, which is completely consis- 
tent with the observed profiles. Virtually all of our best fitting 
models have non-zero value for the star formation lower den- 
sity cutoff A, which results in the outskirts of the galaxy ap- 
pearing to be dynamically colder i n the model <j(r) profile, 
consistent with the obse rvations of IWilkinson et all J2004I) . 
As we showed in § 12.51 the total velocity dispersion in our 
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models with A > does not decrease with radius, with the ob- 
served decline of the line-of-sight velocity dispersion being 
due to significant radial anisotropy exhibited by our models 
outside the initial (non-equilibrium) stellar radius x\ (see Fig- 
ure 

As can be seen from Table |2] the best fitting models for 
Draco have a large range of DM halo virial masses ~ 10 8 - 
10 10 M Q . These masses are all significantly larger than con- 
ventional estimates of ~ 2 x 10 7 M Q (lMateoll998l) which are 
based on the assumption that "mass follows light". The total 
range for halo formation redshifts is ~ 1.3— 10, correspond- 
ing to a galaxy age ~ 9 — 13 Gyr, which is consistent with the 
estim ated age of the oldest Draco stars of 10- 14 Gyr (Mateo 
1998). 

It is interesting to note that a comparable number of both 
NFW and Burkert models provided good fits to the observed 
Draco profiles. (Out of 24 "good" models for Dracol and 
Draco2 datasets, there are 16 NFW and 8 Burkert profiles.) 
Even more interesting is the fact that for 16 "good" NFW 
models v c ~ —1.21 ±0.44, whereas for 8 Burkert models 
v c ~ 1.16 ±0.31. It appears that the Draco galaxy is best 
described by models with either underdense NFW or over- 
dense Burkert DM halos, which is consistent with Draco hav- 
ing a DM halo with an intermediate value of its central density 
slope of 7 ~ -0.5. 

The predictions of our models for the initial gas state 
in Draco appear to be sensible, with the temperature of ~ 
20,000 K (~ 10,000 K if the gas was ionized) and the central 
number density of ~ 10 cm" 3 . The lowest density for star- 
forming gas is ~ 1-2 cm" 3 . After the starburst, the stellar 
cluster in all "good" models expanded by a very modest fac- 
tor of crini.o/oo ~ 1.2. 

Overall, our simple model of single-burst star formation 
from isothermal gas in an extended DM halo appears to pro- 
vide an adequate description of the Draco dSph. 

4.3. Sculptor 

For Sculptor (and als o for Carina and Fo rnax) we used 
star count profiles from W alcher et alJ d2003h . The authors 
use single-color photometry to identify stars belonging to 
these galaxies. As a result, the star count profiles for these 
three dSphs are of much l ower quality than the Draco data of 
lOdenkirchen et all (1200 ll) . who employed a multicolor mask 
filtering technique to improve the sensitivity to Draco stars 
by a factor of ~ 30. One has to keep in mind that sim- 
ple sin gle-color techniques [such as used by IWalcher et al.l 
(2003)] can produce very unreliable results at surface bright- 
ness levels below the foreground/background contamination 
level. For example, according to th e single-color photome- 
try of llrwin & H atzidimit rioul(ll995l) . Draco has a flat surface 
brightness profile beyond 20 arcmin from the galactic cen- 
ter (which was interpreted as a halo of unbound stars). In 
contrast, much higher se nsitivity multicolor observations of 
lOdenkirchen et all(l2001l) failed to find any evidence for such 
a halo. 

In all the three galaxies (Sculptor, Carina, and Fornax) 
some of the radial bin s in their outskirts h ave negative sur- 
face brightness values (Walcher et al. 2003; this is an artifact 
of background subtraction procedure). In the case of Sculp- 
tor, even the upper one-sigma errorbar is negative for the bin 
at r ~ 2 kpc. As a result, we could not show this bin value in 
Figure^. Our \ 2 fitting is not affected by the negative values 
of surface brightness, as we perform the cross-correlation in 
linear space. 



With all the above caveats in mind, our model appears to do 
a reasonably good job in describing the properties of Sculp- 
tor. Formally, we have only one "good" model for this galaxy, 
with x' 2 = 0.94 (see Table [2). On the other hand, we have a 
large number of models with \ 2 slightly worse than for the 
best fitting King model: 3 models with x' 2 < 1-1, ar, d H 
models with x' 2 < 1-5. As will be seen in § 14.51 this is very 
different from the case of Fornax, where the best model has 
x' 2 ~ 1.1, and the second-to-best has x' 2 ~ 1-9. 

Analysis of Figure^ shows that our best fitting model im- 
proves upon theoretical King model by having a less steep 
outer surface bright ness profile, which is the same situation 
as with Draco (see § 14.21 . One cannot have a much better fit 
with a smooth function due to the fact that the data are very 
noisy in the outer parts of the galaxy, with the star counts of 
-0.31 ±0.02 arcmin" 2 at r ~ 2.1 kpc and 0.09 ±0.02 in the 
next radial bin at r ~ 2.6 kpc. 

Properties of three models with x' 2 < 1 ■ 1 are summarized 
in Table |3 With such small statistics it is hard to draw any 
conclusions on the range of model parameters, but it appears 
that all the main parameters for Sculptor are comparable to 
those for Draco. Most importantly, the virial mass for the 
DM halo is required to be > 3 x 10 8 Mb, which is much 
larger than the 6.4 x 10 6 M Q estimate of Mateo ( 1998) based 
on the T = constant assumption. (Interestingly, the two H I 
clo uds observed in the vicinity o f Sculpto r whic h were argued 
by Bouchard. Carignan. & Mashchenko (2003) to be a part 
of ISM of this galaxy, are located in projection well within 
the inferred large tidal radius of ~ 7.5 kpc, corresponding to 
~ 5.5 angular degrees at the distance D = 19 kpc). The no- 
table difference from the Draco case is a factor of two lower 
temperature, which is of order of 10,000 K for neutral gas (see 
Table|2j. 

4.4. Carina 

Similarly to Sculptor, the surface brightnessprofile for Ca- 
rina is very noisy at £ < O.OlEo (see Figures^ and[7t), with 
S being negative at r ~ 1.3 kpc. We found 6 "good" models 
for Carina, with x' 2 = 0.84 for the best fitting model (TableEJ. 
As in the cases of Draco and Sculptor, the surface brightness 
profiles of our Carina models are identical to the best fitting 
theoretical King model in the inner part of the galaxy, and 
become less steep in its outskirts. 

Unlike Draco and Sculptor, all "good" Carina models have 
Burkert DM density profile. All virial masses are larg er than 
~3x 10 s M Q , and are typically around 10 9 M s . | Mateo 
( 1998) gives 1.3 x 10 7 M Q .] Such masses are large enough 
to keep fully photoionized ISM gravitationally bound, which 
is a prerequisite for the radia tion harassment picture proposed 
by Mashch enko et alJ ll2004l) to explain the complex star for- 
mation history of this dwarf (Mateo 1998). All 6 models have 
a halo formation redshift of ~ 5.7 (lookback time 12.6 Gyr), 
which is consistent with the ~ 1 0— 14 Gyr age of the oldest 
stars in this galaxy (Mate"ol [T998h . The rest of the model pa- 
rameters are comparable to the Draco and Sculptor cases. 

4.5. Fornax 

Fornax is the only galaxy in our sample with not a single 
model which woul d sim ultaneously meet all the three selec- 
tion criteria from § 13.11 In Table [2] we list parameters for the 
only model which has a comparable x 2 to the best fitting King 
model: x' 2 = 109. The second-to-best model has x' 2 = 1-9. 

One could try to argue that the failure to find a model with 
x' 2 < 1 is due to the fact that the King model already provides 
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FIG. 6. — Best fitting models for Sculptor, (a) Surface density profiles. The observed profile is shown as solid circles with errorbars. The solid line corresponds 
to the best fitting model (with the lowest x' 2 )- The dashed line shows the best fitting theoretical King model profile, (b) Line-of-sight velocity dispersion profiles 
for all 3 "good" models from Tablel2l The thick line corresponds to the model with the lowest x' 2 - 




FIG. 7. — Best fitting models for Carina, (a) Surface density profiles. The observed profile is shown as solid circles with errorbars. The solid line corresponds 
to the best fitting model (with the lowest x' 2 )- The dashed line shows the best fitting theoretical King model profile, (b) Line-of-sight velocity dispersion profiles 
for all 6 "good" models from Tablel2l The thick line corresponds to the model with the lowest x' 2 - 



an almost perfect fit to the data (see Figure^). It is conceiv- 
able that with unlimited computing resources we would be 
able to fine-tune our initial parameters to find a model which 
would fit the observed surface brightness profile slightly bet- 
ter than the King model. However, King-like profiles are not 
natural for our models (which typically have less steep outer 
density profiles), and it would be very unlikely that it is a 
mere accident that Fornax happened to have an almost perfect 
King profile. The tidal stripping process, on the other hand, 
drives naturally an isothermal stellar system toward a King- 
like state, which makes it a more attractive explanation for 
the Fornax observed properties. 

The inferred parameters for our best fitting model (Ta- 
ble |3 make even more questionable the applicability of our 



model to Fornax. In particular, the inferred central gas den- 
sity pgfl ~ 0.9 cm" 3 is much lower than for the other three 
dSphs. Furthermore, all models with x' 2 iS 3 have A = 0. It 
seems unphysical for star formation to take place with almost 
100% efficiency (see Table|2j at densities <C 1 cm" 3 . 

All of the above leads us to conclude that our model (at least 
in its present, simplest form) is not applicable to the Fornax 
dSph. It remains to be seen if some other non-tidal mecha- 
nism could produce naturally a King-like surface brightness 
profile in a dSph. For now, Fornax being a tidally limited sys- 
tem seems to be the most likely explanation. It is interesting to 
note that even in the tidally limited scenario, Fornax appears 
to have a quite large mass, with the "mass follows light" es- 
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r, kpc r, kpc 

FIG. 8. — Best fitting model for Fornax, (a) Surface density profiles. The observed profile is shown as solid circles with errorbars. The solid line corresponds 
to the best fitting model (with the lowest x' 2 )- The dashed line shows the best fitting theoretical King model profile, (b) Line-of-sight velocity dispersion profile 
for the model with the lowest x' 2 - 



timate of 6.8 x 10 7 M Q ( Mateol [T998T) . and potentially even 
larger value if Fornax is more DM dominated in its outskirts. 

5. DISCUSSION 

The observational data on three galaxies (Draco, Sculptor, 
and Carina) out of four dSphs studied in this paper appear 
to be consistent with our simple model of star formation in 
extended DM halos, and as a consequence are consistent with 
them having very massive DM halos. The galaxy with the best 
quality star count profiles (Draco) is also the one which makes 
the most convinc ing case for having an exte nded halo — both 
in the analysis of Odenki rchen et alJ d2001h who showed that 
the outer stellar isophotes for Draco are not tidally distorted, 
and in the analysis carried out in this paper. 

The data for Fornax appear not to be consistent with the 
predictions of our model. The surface brightness profile for 
this galaxy is almost perfectly described by a theoretical King 
model, which was derived for an isothermal system exposed 
to external tidal field. This type of profile does not come nat- 
urally in our model. In addition, the best fitting model for 
Fornax predicts that the star formation in this dwarf should 
have taken place at unphysically low densities of <C 1 cm -3 . 
The results of our simulations are thus suggestive of the cutoff 
in the outer surface brightness profile of Fornax being caused 
by the tidal field of the Milky Way. 

The situation is significantly different with the three other 
dSphs we studied in this paper (Draco, Sculptor, and Ca- 
rina), which appear to have extended DM halos. This dif- 
ference can be explained if (1) the orbit of Fornax has a sig- 
nificantly smaller pericenter distance than the orbits of the 
other three galaxies [but this seems to be inconsistent with the 
available proper motion measurements (Dinescu et al. 2004), 
which imply that Fornax is on a low eccentricity orbit and is 
currently near its pericenter], and/or (2) the average DM den- 
sity within the stellar extent of Fornax (radius of ~ 4 kpc) is 
significantly smaller than the average DM density within the 
same radius of the three other dSphs. (The latter would be the 
case if Fornax was formed more recently than the three other 



galaxies.) 

The conclusion we reach here is quite contrary to the con- 
ventional argument that the presence of "extratidal" stars 
(those which populate the outer halo of the dwarf galaxy be- 
yond the tidal radius of the best-fitting King model) around 
Milky Way satellites is a good evidence for the dwarf galaxy 
to have undergone a tidal disruption process, with the ob- 
served "extratidal" stars being part of the halo of tidally 
stripped stars. In our model, star formation in extended DM 
halos naturally produces surface brightness profiles which are 
identical to King profiles out to a few core radii, and are less 
steep in the outer parts of the galaxy. The example of Fornax 
suggests that actually the absence of "extratidal" halo, with 
the galactic surface brightness profile being described almost 
perfectly by theoretical King profile, could be an evidence for 
the galaxy to be observed all the way to its tidal radius. In 
such tidally limited systems a small number of stars do have 
to populate an extratidal halo, forming leading and trailing 
tidal tails, but for stable satellite galaxies the surface bright- 
ness of such tails can be very low (especially if dSphs have 
more DM in their outskirts than at the center, which appears 
to be the case for Draco), potentially below the sensitivity of 
existing observations. The systems where tidal tails have been 
reliably detected (such as Sagittarius dSph and Palomar 5) are 
at the last stage of being tidally disrupted, which cannot be the 
case for most Galactic dSphs. We argue that to rule out our 
alternative, non-tidal explanation for the presence of "extrati- 
dal" stars around some Galactic dSphs, unambiguous obser- 
vational evidence for these stars forming leading and trailing 
tidal tails has to be obtained. 

Recent observations suggest that in some Galactic dSphs 
(Draco, Ursa Minor, and Sextans) the very outskirts of the 
galaxies have lower values of line-of-sigh t velocity disper- 
sion than the bulk of the galaxy ( Wilki nson et alJ l2004t 
Klevna et all 120041) . In our models, this is the case for all 
models with sufficiently large lower density cutoffs for star 
formation: A > 0.1 (see Figures^-HJ)). For such models, we 
observe the line-of-sight velocity dispersion to decline by a 
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factor of two in the outer parts of the galaxy, which is caused 
by strong radial anisotropy of stars which were born within 
the initial stellar cluster radius r\ and then later populated the 
r > r\ parts of the galaxy. 

The Fornax (Mat eoll 19971). Dra co (IWilkinson et alJ [2004). 
and Sextans ( Klevna et afl 120041) observations suggest that 
these galaxies are kinematically colder at the center than in 
the intermediate regions. Such behavior cannot be explained 
in terms of our single star burst model. Many dSphs how- 
ever have complex ancient (> 10 Gyr old) star formation his- 
tories, with slightly younger and/or mor e metal rich stars be- 
ing more concentra ted toward the center ( Harbeck et al. 2001 ; 
iTolstov et alJl200l . A natural extension of our model then 
would be to assume that a fraction of gas, which was some- 
what metal-enriched and expelled to large radii during the first 
star burst, was later reaccreted, leading to the second episode 
of star formation. The reaccretion of the gas and the sec- 
ond star burst can be delayed by up to a few Gyrs by su- 
pernovae type la occurring in th e original stellar population 
( Burkert & Ruiz-Lapuente 1997). Let us consider the sim- 
plified case of negligible self-gravity of gas [m' g <C M'(x) in 
eq. 0], with the most of the gas being within the scaling ra- 
dius r s of the DM halo [so that M'(x) oc x a , where a = 2 for 
NFW halos, and 3 for Burkert halos]. With the above simpli- 
fying assumptions, equation (|6j can be rewritten as 



dx 



The solution for this equation is 



p'(x) = p' exp 



Bx' 



a-l 



(21) 



(22) 



where B is some positive constant. For the case of a > 1, the 
total gas mass is 



m' or n' r' 6 ^ a ~ l) 

m g oc p gfi c g 



(23) 

We assume, that a star burst in a slowly cooling and contract- 
ing gas cloud will take place when the central gas density p' 
reaches a certain critical value, which is the same for the orig- 
inal star burst (the total gas mass m' l ) and the subsequent star 
burst (the total gas mass m' 2 ). Using equation d23i . we can 
then relate the gas parameters for the two star bursts through 



4,1 




(a-l)/6 



(24) 



Here c' g l and c' g 2 are the sound speed in the gas in the first 
and second star bursts. The exponent in equation J24i is posi- 
tive for both Burkert and NFW halos (the exponent values are 
1/3 and 1/6, respectively). Due to gas consumption by star 
formation, and potentially some additional gas losses into in- 
tergalactic space, the total gas mass in the second star burst is 
smaller than in the original star burst: m' 2 < m' ; . From equa- 
tion i24\ . this results in the sound speed (and hence the tem- 
perature) of the gas in the second star burst being lower than 
in the first one: c' g2 < c' g l . Both lower temperature and lower 
total gas mass will also result in the gas cloud being more 
compact in the second star burst. As a result, a kinematically 
colder stellar core, slightly younger and slightly more metal 
rich than the bulk of the galaxy, could be formed. A multiple 
star bursts scenario can explain the radial population gradients 



observed in many Local Group dSphs (Har beck et al.l t2001 ). 
and possibly the kink observed in the surface brightne ss pro- 
file o f Draco at the radius of ~ 25 arcmin (Wilkinson et al. 
2004). [N ote that no kink is observed in the Draco profile of 
lodenkirchen et all d200 

For the three galaxies, which appear to be consistent with 
our model (Draco, Sculptor, and Carina), the values of the 
inferred model parameters presented in Table [2] do not differ 
much from one galaxy to another. The "good" models for 
these galaxies have a shallow inner DM density profile (with 
the slope 7 ~ -0.5 ... 0), which is consistent with the observa- 
tions of dwarf spiral galaxies (Burkert 1995), dwarf irregular 
galaxies (ICote. Carignan. & Freema nll2000l). an d low surface 
brightness galaxies (e.g. Marchesi niet al.ll2002l) . The typical 
virial masses for the DM halos are of order of ~ 10 9 M Q , 
with the total range spanning from 9 x 10 7 -2 x 10 10 M©. 
These masses are much larger than those obtained under 
assumption of constant mass-to-light ratio, and are consis- 
tent with the idea that dSphs represent the largest subhalos 
expected to populat e the Milky Way halo in co smological 
ACDM simulations dStoehr et al1l2002t lHavashi et alJl2003l) . 
allevi ating the "missing satellites" problem of such cosmolo- 
gies dMoore et alll999l) . The inferred formation redshifts for 
these galaxies are ~ 2-6 (lookback time ~ 10-13 Gyr), 
which is consistent w ith the age of the oldest stars in these 
galaxies (Mateol ll998l) . The virial radii for the DM halos at 
the time of their formation range from 1 .4 - 30 kpc, whereas 
the tidal radii at their current distance from the Milky Way 
center are 5-32 kpc (typically ~ 10 kpc, which would corre- 
spond to ~ 6 angular degrees in the sky for these three dSphs). 
Both virial and tidal radii are significantly larger than the ob- 
served extent of these galaxies. 

The star forming gas in the "good" models has a tempera- 
ture of- 10, 000-20, 000K if it is neutral (~ 5000- 10,000 K 
if it is ionized), a central density of — 10 cm" 3 , and a central 
pressure of — 10 5 K cm" 3 . The temperature of the gas is a fac- 
tor of — 2.5 lowe r than the virial temperature of the DM halo 
(defined in § !2.4t . suggesting that the gas accreted by the DM 
halo experienced only modest level of radiative cooling glob- 
ally prior to the starburst. The lower density cutoff for star 
formation is of order of — 1 cm" 3 . The inferred star formation 
efficiency in the star-forming central part of ISM is — 10%. 

We argue that our results are not very sensitive to the model 
assumption that the DM potential is static. During the for- 
mation of a dwarf spheroidal galaxy, gradual collapse of the 
gas cloud due to radiative cooling and gravitational instability 
will lead to the adiabatic contraction of the DM halo in the 
central region of the galaxy, where the enclosed gas mass be- 
comes larger than that of DM. Proper modeling of this process 
requires adding many physical processes (live DM halo, non- 
equilibrium chemistry, radiative heating and cooling) with 
many associated free parameters, which is beyond the scope 
of this paper. Instead, we argue that our main results should 
not be significantly affected by the above effect by noting that 
(1) in our best fitting models, gas does not significantly dom- 
inate DM at the center of the galaxy (with the maximum ratio 
of the enclosed gas mass to that of DM of 10 for the Sculptor 
best fitting model; in case of the Draco 1 dataset, this ratio is 
less than 1/3 at any radius), and (2) the initial adiabatic con- 
traction of the DM halo will be compensated by the adiabatic 
expansion caused by the removal of — 90% of baryons from 
the central part of the galaxy by stellar feedback mechanisms 
after the starburst. 

Another potential pitfall is our assumption that the tidal 
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field of Milky Way played no role in shaping the stellar out- 
skirts of the dSphs we studied. Judging from the data pre- 
sented in Table |5] this assumption appears to be justified. In 
particular, the inferred tidal radii are significantly larger than 
the observed extent of these dwarfs (by a factor of > 2). We 
do not expect many dSph stars to be tidally stripped under 
these circumstances. More accurate analysis of this effect is 
complicated as it would involve high resolution iV-body sim- 
ulations of a two-component (DM + stars) dwarf satellite in 
the potential of the host galaxy with some prescription for dy- 
namical friction to describe the orbital decay, and with many 
poorly known orbital and Milky Way halo parameters. This 
task is beyond the scope of our paper, where our main goal 
is to sample well the initial model parameter space, with the 
only practical way of achieving this being in keeping the num- 
ber of free parameters as small as possible. 

The selective nature of the tidal disruption of satellites in 
the halo of the host galaxy could make our method of find- 
in g the most likely DM halo progenitors for Galactic dSphs 
(§ 13.31 invalid. However, this effect would be important 
only if the DM halos in our best fitting models were rela- 
tively easy to disrupt over the Hubble time in the tidal field 
of Milky Way. In the case of the Draco 1 dataset, the best fit- 
ting model has an NFW DM density profile and the inferred 
tidal radius of ~ 4.4 r, (Table |2ji. Numerical simulations of 
Hava shi et alJ J2003) demonstrated that NFW subhalos orbit- 
ing in the host halo potential stay relatively intact after > 10 
orbits if > 2r s . In this respect, the Dracol model appears to 
be stable. In the three remaining datasets (Draco2, Sculptor, 
and Carina; we exclude Fornax from this analysis), the DM 
halos of the best fitting models have Burkert profiles, with the 
tidal radii r t id > 10r s . Burkert satellites are easier to disrupt 
tidally than NFW satellites, as their bin ding radius is ~ 2.1 
times larger (Mashchenko & Sills 2005b). Nevertheless, the 
large inferred value of the tidal radius (> 10r s ) makes it very 
unlikely that the Burkert DM halos corresponding to our best 
fitting models would be completely destroyed by the tidal field 
of Milky Way in a Hubble time. 

We believe that our model is consistent with the observed 
metalli cities o f Galactic dSph s. As can be seen in Fig. 7 of 
Mateo. 1119981) and Fig. 1 of Ta mura et all fcOOll) . the low- 
luminosity objects we are interested in show a relatively large 
spread in metallicity (with a typical value of [Fe/H]^ —1.8) 



and no statistically significant correlation between metallicity 
and luminosity or (conventionally derived) virial mass. The 
only exception among the Galactic dSphs is Fornax, which 
is both the brightest and the most metal-rich ([Fe/H]^ —1.3, 
Mateo Jl998l) . As we discussed above, Fornax is not well 
described by our model. If we relax the single starburst as- 
sumption, more massive galaxies can become more metal-rich 
over time, as they have deeper gravitational potential wells 
and larger tidal radii, enabling them to reaccrete more metal- 
enriched gas expelled in the previous episode of star forma- 
tion. Our model is consistent with the data if the Galactic 
dSphs were formed from a pre-enriched intergalactic medium, 
with [Fe/H]~-2. 

6. CONCLUSIONS 

Our simple model of stars forming in an extended DM halo 
from isothermal gas, and later dynamically relaxing after ex- 
pelling the leftover gas, is consistent with the observed prop- 
erties of three out of four dSphs studied in this paper. 

The results of our simulations suggest that there is an al- 
ternative, non-tidal explanation for the observed presence of 
halos of "extratidal" stars around some Galactic dSphs. In our 
model, such halos are formed when a freshly formed stellar 
cluster is relaxed dynamically inside the extended DM halo. 

The virial masses of DM halos, inferred from fitting our 
models to the three dSphs (Draco, Sculptor, and Carina), are 
in the right range (~ 10 9 M Q ) to give support to the idea that 
the Galactic dSphs represent the most massive subhalos pre- 
dicted to orbit in the Milky Way halo by cosmological A^-body 
ACDM simulations, alleviating in this way the "missing satel- 
lites" problem. The masses are also large enough to keep a 
fully ionized ISM gravitationa lly bound, which is an es sential 
ingre dient of some mod els jBurkert & Ruiz-Laouentdf 19971 
Mashchenko et al. 2004) proposed to explain the complex star 
formation history of some Galactic dSphs. 

We would like to thank Mark Wilkinson and Carl Walcher 
for providing the star count profiles for the dSphs, and Dean 
McLaughlin for making his set of theoretical King models 
available to us. The A^-body simulations reported in this pa- 
per were carried out at the Canadian Institute for Theoretical 
Astrophysics. 



REFERENCES 



Bode, P., Ostriker, J. P., & Turok, N. 2001, ApJ, 556, 93 
Bouchard, A., Carignan, C, & Mashchenko, S. 2003, AJ, 126, 1295 
Boylan-Kolchin, M., Ma, C, & Quataert, E. 2004, ApJ, 613, L37 
Bryan, G. L. & Norman, M. L. 1998, ApJ, 495, 80 

Bullock, J. S., Kolatt, T. S., Sigad, Y., Somerville, R. S., Kravtsov, A. V., 

Klypin, A. A., Primack, J. R., & Dekel, A. 2001, MNRAS, 321, 559 
Burkert, A. 1995, ApJ, 447, L25 
Burkert, A. & Ruiz-Lapuente, P. 1997, ApJ, 480, 297 
Cote, S., Carignan, C, & Freeman, K. C. 2000, AJ, 120, 3027 
Dinescu, D. I„ Keeney, B. A., Majewski, S. R., & Girard, T. M. 2004, AJ, 
128, 687 

El-Zant, A., Shlosman, I., & Hoffman, Y. 2001, ApJ, 560, 636 
Harbeck, D., et al. 2001, AJ, 122, 3092 

Hayashi, E., Navarro, J. E, Taylor, J. E., Stadel, J., & Quinn, T. 2003, ApJ, 
584, 541 

Haywood, M., Robin, A. C, & Creze, M. 1997, A&A, 320, 428 
Irwin, M. & Hatzidimitriou, D. 1995, MNRAS, 277, 1354 
King, I. R. 1966, AJ, 71, 64 

Kleyna, J., Wilkinson, M. I., Evans, N. W., Gilmore, G, & Frayn, C. 2002, 
MNRAS, 330, 792 

Kleyna, J. T., Wilkinson, M. I., Evans, N. W., & Gilmore, G. F. 2004, 

MNRAS, 354, L66 
Lake, G. 1990, MNRAS, 244, 701 
Lokas, E. L. 2002, MNRAS, 333, 697 

Marchesini, D., D'Onghia, E., Chincarini, G., Firmani, C, Conconi, P., 
Molinari, E., & Zacchei, A. 2002, ApJ, 575, 801 



Mashchenko, S., Carignan, C, & Bouchard, A. 2004, MNRAS, 352, 168 
Mashchenko, S., & Sills, A. 2005a, ApJ, 619, 243 
Mashchenko, S., & Sills, A. 2005b, ApJ, 619, 258 

Mateo, M. 1997, ASP Conf. Ser. 1 16: The Nature of Elliptical Galaxies, 259 
Mateo, M. L. 1998, ARA&A, 36, 435 

Mateo, M., Olszewski, E. W., Vogt, S. S., & Keane, M. J. 1998, AJ, 1 16, 2315 
Moore, B., Ghigna, S., Governato, F, Lake, G., Quinn, T, Stadel, J., & Tozzi, 

P. 1999, ApJ, 524, L19 
Navarro, J. E, Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Odenkirchen, M., et al. 2001, AJ, 122, 2538 
Sheth, R. K. & Tormen, G. 1999, MNRAS, 308, 1 19 
Spergel, D. N., et al. 2003, ApJS, 148, 175 

Springel, V., Yoshida, N., & White, S. D. M. 2001, New Astronomy, 6, 79 
Sternberg, A., McKee, C. E, & Wolfire, M. G. 2002, ApJS, 143, 419 
Stoehr E, White S. D. M., Tormen G., Springel V., 2002, MNRAS, 335, L84 
Sugiyama, N. 1995, ApJS, 100, 281 

Tamura, N., Hirashita, H., & Takeuchi, T. T. 2001, ApJ, 552, LI 13 
Tasitsiomi, A. 2003, International Journal of Modern Physics D, 12, 1157 

(preprint astro-ph/0205464 i 
Tolstoy, E. et al. 20U4, ApJ, bl7, LI 19 
van den Bergh, S. 1999, A&A Rev., 9, 273 
van den Bosch, F. C. 2002, MNRAS, 331, 98 

Walcher, C. J., Fried, J. W., Burkert, A., & Klessen, R. S. 2003, A&A, 406, 
847 

Wilkinson, M. I., Kleyna, J., Evans, N. W., & Gilmore, G. 2002, MNRAS, 
330, 778 



Modeling Star Formation in Dwarf Spheroidals 



TABLE 2 
Results of simulations 



Parameter 


Unit 


Statistics 


Draco 1 


Draco2 


Sculptor" 


Carina 


Fornax" 








11 


13 


3 


6 


1 


x' 2 




best 
mean 
range 


0.67 
0.78±0.11 
0.67... 0.97 


0.58 
0.76 ±0.14 
0.58... 1.00 


0.94 
1.03 ±0.08 
0.94... 1.09 


0.84 
0.93 ±0.06 
0.84... 1.00 


1.09 


lgf* 




best 
mean 
range 


0.70 
0.70 ±0.34 
0.02... 1.31 


1.44 
1.27 ±0.69 
0.09... 2.21 


1.89 
0.79 ±0.95 
0.12... 1.89 


1.91 
1.18±0.74 
0.05 . . . 2.07 


1.30 


7 




best 
mean 
range 


-1.00 
-0.82 ±0.40 
-1.00... 0.00 


0.00 
-0.54 ±0.52 
-1.00... 0.00 


0.00 
-0.67 ±0.58 
-1.00... 0.00 


0.00 

0.00 ±0.00 
0.00... 0.00 


0.00 


lg«£ 




best 
range 


-0.79 
_0 7Q 4- O 04 

-0.87... -0.69 


-0.62 

—0 66-1-0 IS 
-1.10. ..-0.49 


-0.80 

_0 79 4- 06 
-0.80... -0.69 


-0.80 
— 1 OS -4- O 93 
-1.40... -0.80 


-0.80 






best 

mean 
range 


1.30 

0.83 ±0.65 
0.00... 1.79 


0.48 

0.66 ±0.68 
-1.00... 1.51 


1.00 

1.36 ±0.42 
1.00... 1.81 


0.00 

-0.12 ±0.44 
-1.00. ..0.14 


0.00 






best 
mean 
range 


-1.00 
-0.66 ±0.45 
-1.10... 0.00 


-1.29 
-0.95 ±0.25 
-1.29. ..-0.41 


-0.63 
-0.37 ±0.28 
-0.63... -0.08 


-1.00 
-0.83 ±0.41 
-1.00... 0.00 


-0.15 


A 




best 
mean 
range 


0.10 
0.14±0.08 
0.10. ..0.30 


0. 10 
0.15±0.11 
0.00... 0.30 


0.00 
0.07 ± 0.06 
0.00... 0.10 


0.10 

0.08 ±0.12 
0.00 . . . 0.30 


0.00 


igM s 


M 


best 
mean 
range 


9.75 
9.47 ±0.33 
8.88... 10.05 


8.48 
8.86 ±0.70 
7.98... 9.96 


8.59 
8.86 ±0.26 
8.59. ..9.12 


8.59 
9.26 ±0.65 
8.49... 10.26 


9.25 


lgr s 


kpc 


best 
mean 
range 


0.54 
0.41 ±0.27 
-0.10... 0.74 


-0.30 
0.02 ±0.42 
-0.56... 0.59 


-0.13 
0.23 ±0.32 
-0.13... 0.47 


-0.12 
0.14 ±0.23 
-0.12... 0.49 


0.27 


LgMvir 


M Q 


best 
mean 
range 


9.79 
9.49 ±0.34 
8.84... 10.06 


8.45 
8.86 ±0.72 
7.93... 9.95 


8.54 
8.84 ±0.28 
8.54. ..9.11 


8.55 
9.23 ±0.66 
8.46... 10.24 


9.27 






best 
mean 
range 


34=^ 1 1 
0.11. ..0.42 


0.10 
2? ± H 
0.09... 0.40 


0.14 
33±0 17 
0.14... 0.43 


l^iJz % 01 
0.14... 0.17 


0.24 






best 
mean 
range 


-0.94 
-0.90 ±1.05 
-1.68... 1.42 


1.16 
-0.01 ±1.23 
-1.60... 1.70 


0.22 
-1.60 ± 1.58 
-2.52... 0.22 


0.14 
0.30 ±0.34 
-0.22... 0.74 


-0.56 


lg(°"ini,o/o-o) 




best 
mean 
range 


0.09 
0.08 ±0.01 
0.06. ..0.10 


0.11 
0.10 ±0.02 
0.07... 0.14 


0.08 
0.09 ± 0.02 
0.08... 0.12 


0.07 
0.05 ± 0.02 
0.03... 0.08 


0.02 


Igf 


K 


best 
mean 
range 


4.43 
4.28 ±0.09 
4.12. ..4.43 


4.35 
4.32 ±0.10 
4. 11... 4.48 


3 93 
3.98 ±0.09 
3.93... 4.08 


3.93 
3.84 ±0.08 
3.76... 3.93 


4.18 


IgPg.O 


cm~ J 


best 
mean 
range 


0.87 
50 ±0 50 
-0.23... 1.04 


1.29 
86 ± 34 
0.10... 1.29 


1.41 
94 ± 48 
0.45... 1.41 


0.38 
16 ±0 46 
-0.77... 0.38 


-0.16 


lgP 


K cnT J 


best 
mean 
range 


5.30 
4.77 ±0.57 
3.89. ..5.41 


5.63 
5.18±0.36 
4.48... 5. 69 


5.35 
4.92 ±0.49 
4.39... 5.35 


4.31 
3.99 ±0.48 
3.03. ..4.31 


4.03 


Pa 


cm -3 


best 
mean 
range 


0.74 
0.66 ±0.54 
0.06... 1.72 


1.94 
1.16±0.90 
0.00... 2.68 


0.00 
0.39 ±0.46 
0.00... 0.90 


0.24 
0.08 ±0.11 
0.00 . . . 0.24 


0.00 


Ig'Vir 


kpc 


best 
mean 
range 


1.31 
1.16±0.31 
0.54... 1.47 


0.35 
0.73 ±0.47 
0.14... 1.36 


0.50 
0.92 ± 0.37 
0.50. ..1.17 


0.51 
0.78 ±0.23 
0.51. ..1.16 


0.99 


Ig^tid 


kpc 


best 
mean 
range 


1.18 
1.10 ±0.07 
1.00... 1.24 


0.88 
0.95 ±0.19 
0.71... 1.35 


0.88 
0.87 ±0.04 
0.83... 0.90 


0.96 
1.18±0.22 
0.91... 1.50 


1.25 



NOTE. — Here N g is the number of "good" models (satisfying all three selection creteria from i5 13 .11 . x' 1 — xVxKina' ^' — ^V^min, 
a = 1/(1 +z), p s .a is the central number density of gas, P = p g .oT is the central gas pressure, p\ = p s fi\, r v ; r is the halo virial radius at the 
formation epoch, and r,id is the current tidal radius (see definition in text). 



"Models with x' 2 < 1.1. 
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